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Abstract 



We compute power corrections to hadronic event shapes in e + e~ annihilation, 
assuming an infrared regular behaviour of the effective coupling a s . With the 
integral of a s over the infrared region as the only non-perturbative parameter, 
also measured in heavy quark physics, we can account for the empirical features 
of l/Q corrections to the mean values of various event shapes. 
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1. Introduction 



Infrared-safe shape measures for hadronic final states in e + e~ annihilation would appear in 
principle to be an ideal testing ground for perturbative QCD. Such quantities are asymp- 
totically insensitive to long-distance non-perturbative physics and can thus be computed 
order-by-order in perturbation theory. The large momentum scale Q ~ Mz available at 
existing e + e~ machines implies a small value of the running coupling a s (Q) and so the 
perturbation series should be relatively well-behaved. Non-perturbative effects should be 
suppressed by inverse powers of Q. Thus hadronic event shapes should be capable of 
providing a high-precision measurement of a s . 

Unfortunately the hoped-for precision has not yet been achieved, partly because 
calculations of event shapes are still lacking, but also because non-perturbative effects turn 
out to be significant even at Q ~ M z . This is because they are in fact suppressed by 
only a single inverse power of Q in many cases. Bearing in mind that a s (Mz) ~ 0.12 
and the non-perturbative scale is 0(1 GeV), we see that the power correction may easily 
be comparable with the 0(a^) next-to- leading term at present energies. Consequently it 
becomes essential to achieve some understanding of power corrections before embarking on 
any 0(a\) calculations of event shapes. 

In the present paper we adopt the approach, advocated in Refs. [|I|,0 , of trying to deduce 
as much as possible about power corrections from perturbation theory. In particular we 
explore the consequences of assuming that a s , defined in some appropriate way, does not 
grow indefinitely at low scales but instead has an infrared-regular effective form. Then 
various moments of a s , integrated over the infrared region, play the role of non-perturbative 
parameters which determine the form and magnitude of power corrections. Since a s is 
supposed to be universal, we obtain relations between the power corrections to various 
observables. 

Our approach is related to that of Korchemsky and Sterman || , and also to several other 
recent papers that discuss power corrections and the ambiguities of perturbation theory in 
terms of infrared renormalons 0, in the context of the Drell-Yan process ||, event shapes 
0, deep inelastic scattering J/J, heavy quark effective theory M or quark confinement M. 
From our viewpoint, infrared renormalons arise from the divergence of the perturbative ex- 
pression for a s at low scales, and the ambiguities associated with different ways of avoiding 
the renormalon poles in the Borel transform plane are resolved by specifying the infrared 
behaviour of a s . This approach implies relationships between the contributions of a given 
renormalon to different processes. 

The quantitative results we obtain look very good in the case of the mean value of 



the thrust parameter pi . The required value of the relevant moment of a s is consistent 



with that obtained from a similar approach to heavy quark fragmentation [TJ . For the other 



shape variables discussed here, the mean value of the C-parameter [nj and the longitudinal 
cross section |12| , a comparison with LEP data is encouraging, but detailed tests must await 
the re-analysis of lower-energy data to establish the energy dependence of the leading power 
correction. 
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2. Calculations 



Consider a quantity of the form 

F=[ Q dkf(k) (1) 

where f(k) behaves like a s (k) k p at k <C Q, say 

f(k) ~ a F a s (A;) P/Q p+1 (A; < Q) (2) 

where we have included the appropriate Q dependence assuming F is dimensionless. Sup- 
pose that F has the perturbative expansion 

F^ t = F 1 a s + F 2 a 2 s + --- . (3) 

More precisely, if the coefficients F n are computed in the MS renormalization scheme at 
scale Q, then in terms of the MS coupling at scale /x R we have 



F^ = F, asOO + [f 2 + ^ In a s 2 (/. R ) + ■ ■ 



(4) 



where /3 = (HCU — 2Nj)/3, with = 3, for Afy active flavours. 

In part of the integration region of Eq. ([[]) the perturbative expression for a s (fc) is not 
appropriate. We may however choose an infrared matching scale fa such that A <C fa <C Q 
and assume that perturbation theory is valid above that scale. We can then introduce a 
non-perturbative parameter a p (fa) to represent the portion of the integral below fa: 

/ dkaJk )¥> = — • (5) 

Jo p + 1 

Before adding this contribution to F pert , we have to subtract the perturbative value of this 
integral, which is, to second order, 



«s(/i R ) + ^ ( In — + — ~~r) «s0 
2tt \ ^ p + 1 J 



R 



(6) 



p + 1 

As a refinement, and for consistency with Ref. 0, we shall assume that the parameter 
atp(fa) refers not to the coupling in the MS scheme but rather to the scheme proposed in 
Ref. fl3 |, which is expected to be more physical in the region under consideration. Thus 



a s in Eq. @ (only) is to be interpreted as af 1 where in terms of the MS coupling 

of = "s + (7) 

with 

K= (*l-*l\c A - b -N l . (8) 
\18 6 j A 9 S K ' 

The only effect on Eq. (|6|) is that the term ln(/i R //i : ) becomes ln(// R / fa) + K/j3 . The full 
expression for F is thus 
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where F pcTt is as given in Eq. (£|) while 



p+l 



(10) 



The dependence of F pert on the renormalization scale /x R is one order higher in a s than that 
calculated, i.e. third-order in this case. Similarly, the dependence of the power correction 
F pow on both // R and the infrared matching scale /x : is third-order, provided /x : is sufficiently 
large for «s(/ii) to have reached its perturbative form. Of course, the value obtained for 
a p (/ii) depends on but this is mostly compensated by the other /independent term. 

The value of the power p and the coefficient ap may be found from the infrared cutoff 
dependence of the lowest-order perturbative result. In this connection, it is crucial that the 
appropriate argument of a s for soft and/or collinear gluon emission is the gluon transverse 



momentum k i [141 . Thus the cutoff should be a k -cutoff. 



Consider for example the mean value of the thrust T. The contribution to this quantity 
from the region k±_ < /x r is 

6(T) = -— dx 1 dx 2 a s (k ± )- 1 2 r min{(l - a?i), (1 - x 2 )} (11) 

2n Jk ± <m (1 — Xi){l — x 2 ) 

where Cp = 4/3. Setting 1 — xi j2 = yi j2 an d considering the region yi < y 2 <C 1, we have 
k± = \Fy~\!yiQ and hence 

r ,ms C F m/Q , m dk ± , 

8{T) = -4-^ / tZj/! / -±a s (k L ) 

7T JO ^s/i Q k±_ 

r d k ± a s (k ± ) = -4^^« (/ii). (12) 

Thus in this case p = and we obtain a 1/Q correction, with a coefficient in Eq. ( |10D of 
of = ax where 

a T = -4— = -1.70 . (13) 

7T 

As shown by the solid curve in Fig. 1, an excellent fit to the data on (T) over the range 
14 < Q < 92 GeV can be obtained using the perturbative prediction |n| 

(T) pcrt = 1 - 0.335 a s - 1.02 a 2 (14) 

with (J, R = Q and a s (Mz) = 0.117 ± 0.005 [[HJ, plus a power correction of the form QTIj). 
For fii = 2 GeV, the fitted value of the non-perturbative parameter «o is 

r2 GeV 

a (2 GeV) = (2 GeV)" 1 / dkaf(k) = 0.53 ± 0.04 , (15) 

Jo 

with x 2 — 4.5 for 8 degrees of freedom. Allowing both a s (Mz) and do{2 GeV) to be free 
parameters gives 

a 3 (M z ) = 0.120 ±0.004 , a (2 GeV) = 0.52 ± 0.03 , (16) 
with rVd-o.f. = 3.7/7. 
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Figure 1: Mean value of 1 — T, where T is the thrust. 



We also obtain good fits for other values of the arbitrary infrared matching parameter 
fa. At fa = 3 GeV, for example, we find 

a 3 (M z ) = 0.118 ±0.004 , a (3 GeV) = 0.42 ± 0.03 , (17) 

with x 2 /d.o.f. = 4.0/7. The change in a implies that ag ff (2.5 GeV) ~ 0.2 ± 0.1, which 
is reasonable, the perturbative value being around 0.3. The change in the overall power 
correction is small (about 5%), since, as explained above, the /i r dependence mostly cancels 
in Eq. flll]). 

Using the value (p~5| ) of «o, obtained by fitting the thrust data, one can now predict 
the power corrections to other event shapes. For the mean value of the C-parameter, for 
example, we find that the coefficient of in Eq. ( |T(]| ) is 

a c = 6C F = 8. (18) 

At present there are only data on (C) at Q = Mz, where Eqs. (|15|) and (|18| ) imply 

(C) pow = 0.057 ± 0.008 . (19) 



The second-order perturbative prediction is fl5 | 

(C) pcTt = 1.375 a s + 3.88 a 2 s = 0.214 ± 0.011 (20) 
for ct s = 0.117 ± 0.005. Thus the full theoretical prediction is 

(C) th = (C) pcrt + (C) pow = 0.271 ± 0.014 , (21) 
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which is consistent with the experimental result 
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(C) cxp = 0.2587 ± 0.0013 ± 0.0018 . (22) 

Note that the power correction represents over 20% of the value of this quantity. 
Finally, for the longitudinal cross section fraction a h /a tot we predict a coefficient 

a L = C F = 1.33 , (23) 

leading to the power correction at Q = Mz 

K/a tot ) pow = 0.010 ± 0.001 . (24) 

The first-order perturbative prediction is a s /ir = 0.037. However, the second-order correc- 
tion is not yet known. The preliminary OPAL measurement |18|] is 



K/^tot)° xp = 0.067 ± 0.008 . (25) 

Clearly the second-order perturbative correction should be significant if there is to be 
satisfactory agreement between theory and experiment. 

The values for q.q obtained above from event shapes are in reasonable agreement with 
those deduced from heavy quark fragmentation spectra. In Ref. [0 the value cxq{2 GeV) ~ 
0.6 was obtained from fits to heavy quark energy losses in e + e~ annihilation. The same 
conclusion follows from an analysis of the quantity = — In (xh), where xh is the energy 
fraction carried by the heavy quark H, using the approach of the present paper. We find 
a quark mass (1/M) correction of the form (pT0|) , with Q replaced by M, /i R ~ M, and 
coefficient a H = Cp/2. The perturbative prediction deduced from Ref. Q is 

+ -}-(K + 5 2 )[a s (M)-a a (Q)]) , (26) 
Po J 

where 5 2 is the (numerically negligible) 2-loop anomalous dimension correction 

= y) + - = -°- 173 • ( 2? ) 

The expression (p6|), which accounts for the a s \n(Q/M), a s and a|ln(Q/M) terms in £ H , 
but neglects a|? terms with no large logarithm, gives £f rt = 0.26 ± 0.02 for 6-quarks at 
Q = Mz- Comparing with the experimental value of 0.36 ± 0.02 deduced from lepton 
spectra this implies that ^ pow = 0.10 ± 0.03 and hence that a (2 GeV) ~ 0.6 ± 0.1. 



The errors are estimated conservatively, taking into account the small scale ~ 5 GeV 
and the lack of a complete C(a|) calculation of ££ ert . 
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3. Conclusions 



Note that the power correction coefficients ax, ac an d &l deduced above using a k± cutoff 
are identical to those obtained in Ref. f|| with a gluon mass cutoff. With a k± cutoff, 
however, these coefficients have a physical interpretation: they measure the contribution 
of the low-scale region in which a s departs significantly from its perturbative form. After 
being used to calculate the coefficients, the cutoff is replaced by an infrared matching 
parameter which represents the scale below which we switch from the perturbative to 
the non-perturbative description of a s . As long as \l x is not too small (larger than about 
1 GeV) the predictions are quite insensitive to its value, indicating that the perturbative 
behaviour has set in at that scale. 

The divergence in the perturbative expression for a s at low scales is responsible for the 
divergence of the perturbation series for quantities like those considered here, giving rise to 
the so-called "renormalon ambiguity" . By assuming an infrared regular form for the effec- 
tive coupling, we resolve this ambiguity, at the price of introducing the non-perturbative 
parameters a p . These parameters are, however, universal, and can be measured experi- 
mentally, like «o in Eq. (15|). 



Combined fits to the non-perturbative parameters a p and the perturbative parameter 
a s , using data on several different event shapes, provide the possibility of understanding 
something new about QCD at low scales and at the same time measuring a s with improved 
precision. This would be useful not only for QCD but also in constraining physics beyond 
the Standard Model. 
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